Perturbing Conformal Turbulence

We consider perturbations of the non-unitary minimal model solutions of two-dimensional conformal turbulence proposed by Polyakov. Demanding the absence of non-integrable singularities in the resulting theories leads to constraints on the dimension of the perturbing operator. We give some general solutions of these constraints, illustrating with examples of specific models. We also examine the effect of such perturbations on the Hopf equation and derive the interesting result that the latter is invariant under a certain class of perturbations, to first order in perturbation theory, examples of which are given in specific cases.Comment: uuencoded ps file, 13 page

A review of urban roughness sublayer turbulence

It is becoming increasingly important that we can understand and model flow processes in urban areas. Applications such as weather forecasting, air quality and sustainable urban development rely on accurate modelling of the interface between an urban surface and the atmosphere above. This review gives an overview of current understanding of turbulence generated by an urban surface up to a few building heights, the layer called the roughness sublayer (RSL). High quality datasets are also identified which can be used in the development of suitable parameterisations of the urban RSL. Datasets derived from physical and numerical modelling, and full-scale observations in urban areas now exist across a range of urban-type morphologies (e.g. street canyons, cubes, idealised and realistic building layouts). Results show that the urban RSL depth falls within 2 – 5 times mean building height and is not easily related to morphology. Systematic perturbations away from uniform layouts (e.g. varying building heights) have a significant impact on RSL structure and depth. Considerable fetch is required to develop an overlying inertial sublayer, where turbulence is more homogeneous, and some authors have suggested that the “patchiness” of urban areas may prevent inertial sublayers from developing at all. Turbulence statistics suggest similarities between vegetation and urban canopies but key differences are emerging. There is no consensus as to suitable scaling variables, e.g. friction velocity above canopy vs. square root of maximum Reynolds stress, mean vs. maximum building height. The review includes a summary of existing modelling practices and highlights research priorities

Detection of sea-breeze events around London using a fuzzy-logic algorithm

We present an algorithm for detecting sea breezes based on fuzzy logic, using changes in variables commonly measured at meteorological stations. The method is applied to one year’s worth of UK Met Office data (2012) measured at several stations around London, UK. Results indicate about a dozen potential events over the year, when matched against corresponding detections at a coastal reference site (Gravesend). In some cases the time lags between corresponding events detected at different stations can be used to characterise the average propagation speed of the sea-breeze front. Advantages and disadvantages of the method are discussed

Dispersion of a passive scalar within and above an urban street network

The transport of a passive scalar from a continuous point-source release in an urban street network is studied using direct numerical simulation (DNS). Dispersion through the network is characterized by evaluating horizontal fluxes of scalar within and above the urban canopy and vertical exchange fluxes through the canopy top. The relative magnitude and balance of these fluxes are used to distinguish three different regions relative to the source location: a near-field region, a transition region and a far-field region. The partitioning of each of these fluxes into mean and turbulent parts is computed. It is shown that within the canopy the horizontal turbulent flux in the street network is small, whereas above the canopy it comprises a significant fraction of the total flux. Vertical fluxes through the canopy top are predominantly turbulent. The mean and turbulent fluxes are respectively parametrized in terms of an advection velocity and a detrainment velocity and the parametrization incorporated into a simple box-network model. The model treats the coupled dispersion problem within and above the street network in a unified way and predictions of mean concentrations compare well with the DNS data. This demonstrates the usefulness of the box-network approach for process studies and interpretation of results from more detailed numerical simulations

Resolution dependence of turbulent structures in convective boundary layer simulations

Large-eddy simulations are performed using the UK Met Office Large Eddy Model to study the effects of resolution on turbulent structures in a convective boundary layer. A standard Smagorinsky subgrid scheme is used. As the grid length is increased the diagnosed height of the boundary layer increases and the horizontally- and temporally-averaged temperature near the surface and in the inversion layer increase. At the highest resolution, quadrant analysis shows that the majority of events in the lower boundary layer are associated with cold descending air, followed by warm ascending air. The largest contribution to the total heat flux is made by warm ascending air, with associated strong thermals. At lower resolutions, the contribution to the heat flux from cold descending air is increased and that from cold ascending air is reduced in the lower boundary layer; around the inversion layer, however, the contribution from cold ascending air is increased. Calculations of the heating rate show that the differences in cold ascending air are responsible for the warm bias below the boundary layer top in the low resolution simulations. Correlation length and time scales for coherent resolved structures increase with increasing grid coarseness. The results overall suggest that differences in the simulations are due to weaker mixing between thermals and their environment at lower resolutions. Some simple numerical experiments are performed to increase the mixing in the lower-resolution simulations and to investigate backscatter. Such simulations are successful in reducing the contribution of cold ascending air to the heat flux just below the inversion, although the effects in the lower boundary layer are weaker

Flow structure and near-field dispersion in arrays of building-like obstacles

Dispersion in the near-field region of localised releases in urban areas is difficult to predict because of the strong influence of individual buildings. Effects include upstream dispersion, trapping of material into building wakes and enhanced concentration fluctuations. As a result, concentration patterns are highly variable in time and mean profiles in the near field are strongly non-Gaussian. These aspects of near-field dispersion are documented by analysing data from direct numerical simulations in arrays of building-like obstacles and are related to the underlying flow structure. The mean flow structure around the buildings is found to exert a strong influence over the dispersion of material in the near field. Diverging streamlines around buildings enhance lateral dispersion. Entrainment of material into building wakes in the very near field gives rise to secondary sources, which then affect the subsequent dispersion pattern. High levels of concentration fluctuations are also found in this very near field; the fluctuation intensity is of order 2 to 5

Processes controlling atmospheric dispersion through city centres

We develop a process-based model for the dispersion of a passive scalar in the turbulent flow around the buildings of a city centre. The street network model is based on dividing the airspace of the streets and intersections into boxes, within which the turbulence renders the air well mixed. Mean flow advection through the network of street and intersection boxes then mediates further lateral dispersion. At the same time turbulent mixing in the vertical detrains scalar from the streets and intersections into the turbulent boundary layer above the buildings. When the geometry is regular, the street network model has an analytical solution that describes the variation in concentration in a near-field downwind of a single source, where the majority of scalar lies below roof level. The power of the analytical solution is that it demonstrates how the concentration is determined by only three parameters. The plume direction parameter describes the branching of scalar at the street intersections and hence determines the direction of the plume centreline, which may be very different from the above-roof wind direction. The transmission parameter determines the distance travelled before the majority of scalar is detrained into the atmospheric boundary layer above roof level and conventional atmospheric turbulence takes over as the dominant mixing process. Finally, a normalised source strength multiplies this pattern of concentration. This analytical solution converges to a Gaussian plume after a large number of intersections have been traversed, providing theoretical justification for previous studies that have developed empirical fits to Gaussian plume models. The analytical solution is shown to compare well with very high-resolution simulations and with wind tunnel experiments, although re-entrainment of scalar previously detrained into the boundary layer above roofs, which is not accounted for in the analytical solution, is shown to become an important process further downwind from the source

Turbulence characteristics across a range of idealized urban canopy geometries

Good representation of turbulence in urban canopy models is necessary for accurate prediction of momentum and scalar distribution in and above urban canopies. To develop and improve turbulence closure schemes for one-dimensional multi-layer urban canopy models, turbulence characteristics are investigated here by analyzing existing large-eddy simulation and direct numerical simulation data. A range of geometries and flow regimes are analyzed that span packing densities of 0.0625 to 0.44, different building array configurations (cubes and cuboids, aligned and staggered arrays, and variable building height), and different incident wind directions (0 degrees and 45 degrees with regards to the building face). Momentum mixing-length profiles share similar characteristics across the range of geometries, making a first-order momentum mixing-length turbulence closure a promising approach. In vegetation canopies turbulence is dominated by mixing-layer eddies of a scale determined by the canopy top shear length scale. No relationship was found between the depth-averaged momentum mixing length within the canopy and the canopy top shear length scale in the present study. By careful specification of the intrinsic averaging operator in the canopy, an often-overlooked term that accounts for changes in plan area density with height, is included in a first-order momentum mixing-length turbulence closure model. For an array of variable-height buildings, its omission leads to velocity overestimation of up to 17%. Additionally, we observe that the von Karman coefficient varies between 0.20 and 0.51 across simulations, which is the first time such a range of values has been documented. When driving flow is oblique to the building faces, the ratio of dispersive to turbulent momentum flux is larger than unity in the lower half of the canopy, and wake production becomes significant compared to shear production of turbulent momentum flux. It is probable that dispersive momentum fluxes are more significant than previously thought in real urban settings, where the wind direction is almost always oblique